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tin: Abstract 

■ We introduce a family of directed geometric graphs, denoted G\, that depend on two pa- 

rameters A and 6. For < < f and ^ < A < 1, the graph is a strong <-spanner, with 
t — (i-A) cosg - '^^^ out-degree of a node in the graph is at most [27r/ min((?, arccos ^)J . 
Moreover, we show that routing can be achieved locally on G\ . Next, we show that all strong 
t-spanners are also <-spanners of the unit disk graph. Simulations for various values of the pa- 
, rameters A and 6 indicate that for random point sets, the spanning ratio of G^ is better than the 
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proven theoretical bounds. 
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1 Introduction 



A graph G whose vertices are points in the plane and edges are segments weighted by their length 
is a geometric graph. A geometric graph G is a t-spanner (for t > 1) when the weight of the shortest 
^ I path in G between any pair of points a, b does not exceed t ■ \ab\ where \ab\ is the Euclidean distance 

O ' between a and b. Any path from a to 6 in G whose length does not exceed t ■ \ab\ is a t-spanning 

■ path. The smallest constant t having this property is the spanning ratio or stretch factor of the graph. 

A f-spanning path from a to 6 is strong if the length of every edge in the path is at most \ab\. The 
graph G is a strong t-spanner if there is a strong t-spanning path between every pair of vertices. 
^ ■ The spanning properties of various geometric graphs have been studied extensively in the liter- 

^ ! ature (see the book by Narasimhan and Smid |16I] for a comprehensive survey on the topic) . We are 

particularly interested in spanners that are defined by some proximity measure or emptiness crite- 
rion (see for example Bose et al. [1]). Our work was initiated by Chavez et al. [4] who introduced 
a new geometric graph called Half-Space Proximal (HSP). Given a set of points in the plane, HSP is 
defined as follows. There is an edge oriented from a point p to a point q provided there is no point 
r in the set that satisfies the following conditions: 

1. \pr\ < \pq\, 

2. there is an edge from p to r and 

3. g is closer to r than to p. 
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The authors show that this graph has maximum out-degreq^ at most 6. The authors also claim 
that HSP has an upper bound of 27r + 1 on its a stretch factoiH and that this bound is tighld. 
Unfortunately, in both cases, we found statements made in their proofs of both the upper and 
lower bounds to be erroneous or incomplete as we outline in Section[3l However, in reviewing their 
experimental results as well as running some of our own, although their proofs are incomplete, we 
felt that the claimed results might be correct. Our attempts at finding a correct proof to their claims 
was the starting point of this work. Although we have been unable to find a correct proof of their 
claims, we introduce a family of directed geometric graphs that approach HSP asymptotically and 
possess several other interesting characteristics outlined below. 

In this paper, we define a family of graphs. These are directed geometric graphs that depend 
on two parameters A and 6. We show that each graph in this family has bounded out-degree and 
is a strong t-spanner, where both the out-degree and t depend on A and 6. Furthermore, graphs in 
this family admit local routing algorithms that find strong t-spanning paths. Finally, we show that 
all strong t-spanners are also spanners of the unit-disk graph, which are often used to model adhoc 
wireless networks. 

The remainder of this paper is organized as follows. In Section [H we introduce the graph 
and prove its main theoretical properties. In Section [3l we compare the G\ graph to HSP. In 
Section|4l we show that by intersecting the graph with the unit disk graph, we obtain a spanner 
of the unit disk graph. In Section [5l we present some simulation results about the G\ graph. 

2 The family G{ of graphs 

In this section, we define the G\ graph and prove that it is a strong i-spanner of bounded out- 
degree. We first introduce some notation. Let P be a set of points in the plane, < ^ < ^ and 
i<A<l. 

Definition 2.1 The 9-cone{p, r) is the cone of angle 26 with apex p and having the line through p and 
r as bisector. 

Definition 2.2 The A-half-plane(p, r) is the half-plane containing r and having as boundary the line 
perpendicular to pf and intersecting pf at distance ^ \pr\ from p. 

Definition 2.3 The destruction region of r with respect to p, denoted K{p, r), is the intersection of 
the 6-cone{p, r) and the X-half-plane{p, r) (see Figure\l]). 

The directed graph G^^{P) is obtained by the following algorithm. For every point p e P,do the 
following: 

1. Let N{p) be the set P\{p}. 

2. Let r be the point in N{p) which is closest to p. 

3. Add the directed edge {p,r) to 

4. Remove all q G K{p, r) from N{p) (i.e., N{p) ^ N{p) \ K{p, r)). 

5. If N{p) is not empty go to El 

^Theorem 1 in |4] 
^Theorem 2 in |4] 
^Construction in Figure 2 in Hi 
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Figure 1 : The destruction region of r with respect to p. 



The graph computed by this algorithm can ahernatively be defined in the following way: 

Definition 2.4 The directed graph is the graph having P as vertex set and there is an edge 

{p,q) G G^\{P) iff there is no point r e P, such that \pr\ < \pq\, {p,r) is an edge of G^^{P) and 
q G K{p, r), (ties on the distances are broken arbitrarily). Such a point r is said to he a destroyer of 
the edge {p, q). 

2.1 Location of Destroyers 

What prevents the directed pair {p, q) from being an edge in G^? It is the existence of one point 
acting as a destroyer. Given two points p, q, where can a point lie such that it acts as the destroyer 
of the edge {p,q)'^ In this subsection, we describe the region containing the points r such that 
q G K{p,r). This region is denoted K{p,q). In other words, K{p,q) is the description of all the 
locations of possible destroyers of an edge {p, q). 



Proposition 2.5 Let R{p, q, A) be the intersection of the disks Ci centered at p with radius \pq\ and 
C2 centered at c = p + \{q — p) with radius X\pq\. If q ^ K{p, r) and \pr\ < \pq\, then r £ R{p, q, A). 

Proof: If r destroyed {p, q), then \pr\ < \pq\. Therefore, r is in Ci. To complete the proof, we need 
to show that r is in C2. We begin by considering the case when q lies on the line / which is the 
boundary of A-half-plane(p, r) (see Figure [2]). Let si be the midpoint oipf, ti the intersection of I 




Figure 2: The location of a point r destroying the edge {p, q). 



3 











^~ ^ —^=(0,0) 



q p 











^ --^=(0,0) 



A|pg| 
Case 2.2 



Figure 3: Cases for the proof of Theorem |2.8[ 

with pf and c' the intersection of pq with the bisector of pr. Since the triangles Aptiq and Apsic' 
are similar, this implies that 

Iptil 2|pti| 2\pr\ 

Therefore, c' = c, which implies that \cr\ = \cp\ thereby proving that r is on the boundary of C2. 

In the case when q is not on then we have \pc'\ < \pc\ and r lies on a circle centered at c' going 
through p. Therefore, r is contained in C2, which completes the proof. □ 



The following proposition follows directly from the definition of K{p, r). 

Proposition 2.6 Ifq£ K{p,r), then Zqpr < 6. 

Combining Proposition 12.51 and Proposition |2.6[ we get: 

Proposition 2.7 Let K{p, q) be the intersection of R{p, q, A) with the 0-cone{p, 
\pr\ < \pq\, then r £ K{p, q). 



)■ If q £ K{p, r) and 



2.2 The Stretch Factor of 



Theorem 2.8 For < 6* < f and ^ < X < 1, the graph is a strong t-spanner, with t 



(1- A) cose- 
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Proof: Let P be a set of points in the plane, p,q e P and ddp, q) be the length of the shortest path 
from p to g in G^^{P). We show by induction on the rank of the distance \pq\ that ^^(p, q) < t\pq\. 

Base case: If p and q form a closest pair, then the edge {p, q) is in G^^{P) by definition. Therefore, 

dG{p,q) = \pq\ < t\pq\. 

Inductive case: If the edge {p,q) is in G\{P), then dG{p,q) = \pq\ < t\pq\ as required. We now 
address the case when {p, q) is not in G^^{P). By Proposition 12.71 there must be a point r G K{p, q) 
with \pr\ < \pq\ that is destroying {p, q) and such that the edge {p, r) is in G^(P). Since r G K{p, q) 
and \pr\ < \pq\, we have that \rq\ < \pq\. By the inductive hypothesis, we have dcir, q) < t\rq\. 

Let z be the intersection of the boundaries of the disks Ci and C2 defined in Proposition 12. 51 We 
assume, without loss of generality, that c is the origin and that points p, q are on the x-axis with p 
to the left of q as depicted in Figure [3j The remainder of the proof addresses two cases, depending 
on whether or not < (the notation px denotes the x-coordinate of a point p). 

Case 1: < Zx- Let v G K{p,q) be the point with the same x-coordinate as r and having the 
greatest y-coordinate. In other words, v is the highest point in K{p, q) that is strictly above r. We 
have: 

dG{p,q) < \pr\+dG{r,q) 

< \pr\ +t|rq'|(ind. hyp.) 

< \pv\ + t\vq\. 

Now, let a = Zvpq < 9 We express \pv\ and \vq\ as a function of cos a. Consider the triangle A{pvc) 
and note that \vc\ = \pc\ by construction. We have 

\pv\ = 2X\pq\ cos a 

and, from the law of cosines, 

l^^^l^ = Ipt"!^ + Ipgp — 2|pt;| cos a 

= 4A^|p(7p cos^ a + \pq\'^ — 4A|pg|^ cos^ a 
= |M|^(4A^cos^a - 4Acos^a + 1) 

which implies that: 

dciPiQ) < 2A|pg| cos a + t\pq\ \l 4A^ cos^ a — 4A cos^ a + 1 
= |p(7|(2Acos a + t\/ cos^ a — 4Acos^ a + 1). 

Therefore, we have to show that: 

t > 2A cos a + t\J 4A^ cos^ a — 4A cos^ a + 1, 

which can be rewritten as 

^ ^ 2A cos a 

1 — \/ 4A2 cos^ a — 4A cos^ a + 1 
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Since a <9 <tt/2 implies cos 6 < cos a, by straightforward algebraic manipulation we have that 



1 2A cos a 

> 



(1- A) cos a 1 - V4A2 cos^ a - 4A cos^ a + l 

Case 2: rx > Zx. Let p = Zzpq. We first compute the value of cos /?. From the definition of Ci and 
C2, we have 

and 

[Zx - Px) + Zy = \pq\ . 
Therefore, since px = —X\pq\,^e have Zx = ^^^^^2A^^ which implies 

\\pq\ + Zx , , 1-2A^ 1 

cos/3 = ; ; = A H — = — . 

^ \pq\ 2A 2A 

We need to consider two subcases, depending on whether or not (3 <9. 

Case 2.\: (3 <9.ln this case, we have: 

dG{p,q) < \pr\+dG{r,q) 

< \pr\ +t|rg|(ind. hyp.) 

< \pz\ + t\zq\ 

= \pq\ + t\zq\- 

By the law of cosines, we have 

\zqf = |pg|2(2-i) 
which implies 

dGip,q) < \pq\+t\pq\\2 ^ 



A' 



Therefore, we have to show that 
1 



t > 



2-i 



Since P < 9,'we have cos /3 > cos 9, and thus 
1 

t = 
> 



(1- A) cos9 
1 

(1 - A) cos (3 
1 

(1-A)(1/2A) 

2A 
1-A 



Figure 4: The graph has bounded out-degree. 

where the last inequahty holds because it is equivalent to (1 — A)^ > 0. 

Case 2.2: P > 9. By the law of cosines we have 

dciPiQ) < ImI + V2 — 2 cos 6 

which means that we have to show that 
1 

t > , 

1 - \/2 - 2 cos 

But - — 7=i= — = > 71 — 4 — Q since 13 > 6 implies cos 6 > cos 3 = 4\- This completes the last case of 

1-^2-2 008 6* — (1- A) cose I- t- I" 2\ 

the induction step. Note that the resulting t-spanning paths found in this inductive proof are strong 
since both \pr\ and \rq\ are shorter than \pq\. □ 

The above proof provides a simple local routing algorithm. A routing algorithm is considered 
local provided that the only information used to make a decision is the 1 -neighborhood of the 
current node as well as the location of the destination (see [3] for a detailed description of the 
model). The routing algorithm proceeds as follows. To find a path from p to q, if the edge (p, q) is 
in G^-y{P), then take the edge. If the edge (p, q) is not in G^-y{P), then take an edge (p, r) where r is 
a destroyer of the edge (p, q). Recall that r is a destroyer of the edge (p, g) if r G K{jp^ q). This can 
be computed solely with the positions of q and r. Therefore, determining which of the neighbors 
of p in G^^{P) destroyed the edge (p, q) is a local computation. 

2.3 is of Bounded Out-Degree 

Proposition 2.9 The out-degree of a node in the graph is at most [27r/ m\n{6, arccos ^)J. 

Proof: Let {p, r) and (p, s) be two edges of the G\ graph. Without loss of generality, \pr\ < \ps\. Let 
/ be the line perpendicular to pr through p+ -^{r — p). Then either Zspr > 6* or s lies on the same 
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side of / as p. In the latter case, the angle Zspr is at least arccos 2^ (see Figure |4l). The angle Zspr 
is then at least min(6', arccos which means that p has at most [27r/ min(6', arccos ^)J outgoing 
edges. □ 

Corollary 2.10 If 6 > vr/S and A > 2cos(27r/7) ' ^^^^ out-degree of a node in the graph is at 
most six. 

3 Half-Space Proximal 

In this section, we outline the inconsistencies within statements of the proof of the upper and lower 
bounds of HSP given in Chavez et al. [1^. 




Figure 5 : Counter-example to the proof of Theorem 2 of Hi. 

In the proof of the upper bound (Theorem 2 of Chavez et al. [4]), claim 4 states that all ver- 
tices no, Ul, U2, . . . , Ufc are either in clockwise or anticlockwise order around v. The claim is that this 
situation exists when {u, v) is not an edge of HSP and no neighbor of u is adjacent to v. However, 
as stated, this is not true. A counter-example to this claim is shown in Figure [5l There is a unique 
path from u to v, namely uuiU2V, but this path is neither clockwise nor counter-clockwise around v. 
We believe that this situation may exist in the worst case. However, a characterization of the worst 
case situation must be given and it must be proven that the worst case situation has the claimed 
property. 




Figure 6: The illustration of the lower bound on the spanning ratio of f^. 

For the lower bound, the authors also claim that the stretch factor of HSP can be arbitrarily 
close to 27r + 1. However, the proof they provide to support that claim is a construction depicted in 
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Figure [6] (reproduced from iQ] ) ■ The claim is that the path from utov can have length arbitrarily 
close to (27r + l)|nv|. Although this may be true for the path that they highlight. This path is 
not the only path from n to ?; in HSP. The authors neglected the presence of the edge (uiUk) in 
their construction, which provides a shortcut that makes the distance between u and v much less 
than 2\uv\. 

One of the main reasons we believe the claims made in Chavez et al. may be true is that in 
the simulations, all the graphs have small stretch factor. In fact, the stretch factor seems to be even 
smaller than 27r + 1. However, at this point, no proof that HSP is a constant spanner is known. We 
provide a lower bound of 3 — e on the stretch factor of HSP as depicted in the construction below. 



1 - 2S 




Figure 7: Example of a 6 nodes HSP with a stretch factor of 3 — e. The solid edges are in HSP. 



Proposition 3.1 The HSP graph has stretch factor at least 3 — e. 

Proof: Consider the set of 6 point as in Figure [71 put S = e/6. The length of the path between p 
and q via a and b is equal to the length of the path between p and q via c and d. The length of both 
of these paths is 3 — 66. Since the shortest path between p and q in the HSP graph is one of the 
above paths, the stretch factor is 3 — 65 = 3 — e. □ 



4 Unit Disk Graph Spanners 

In Section [21 we showed that the graph of a set of points in the plane is a strong t-spanner 
of the complete graph of these points, for a constant t = cose ■ show in this section that 
strong t-spanners are also spanners of the unit disk graph. That is, the length of the shortest path 
between a pair of points in the unit disk graph is not more than t times the length of the shortest 
path in the graph resulting from the intersection of a strong t-spanner the a unit disk graph. Before 
proceeding, we need to introduce some notation. 

For simplicity of exposition, we will assume that given a set P of points in the plane, no two 
pairs of points are at equal distance from each other. The complete geometric graph defined on a set 
P of points, denoted C{P), is the graph whose vertex set is P and whose edge set is P x P, with 
each edge having its weight equal to the Euclidean distance between its vertices. Let ei, . . . , e^n-j 

be the edges of C{P) sorted according to their lengths Li, . . . , L^,i^ . For i = 1 ... (2), we denote by 
Ci{P) the geometric graph consisting of all edges whose length is no more than Lj. In general, for 
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any graph G whose vertex set is V, we define as G n Ci{V). Let UDG(P) be the unit disk graph 
of P, which is the graph whose vertex set is P and with edges between pairs of vertices whose 
distance is not more than one. Note that UDG(P) = Ci{P) for some i. 

We now show the relationship between strong t-spanners and unit disk graphs. 

Proposition 4.1 If S is a strong t-spanner of C{P), then for all i = I . . .{^ and all j = 1. . .i, Si 
contains a t-spanning path linking the vertices ofej. 

Proof: Let p and q be the vertices of ej. Consider a strong t-spanner path in S between p and q. 
Each edge on this path has length at most \pq\ = Lj < Lj. Therefore, each edge is in 5j. □ 



Proposition 4.2 If S is a strong t-spanner of C{P), then for all i = I. . . (^), Si is a t-spanner of 

Ci{P). 

Proof: Let a and h be any two points such that dQ.(^p^{a, b) is finite. We need to show that in Si 
there exists a path between a and b whose length is at most t ■ dc^[p){a, b). Let a = pi,p2, ■ ■ ■ ,Pk = b 
be a shortest path in Ci{P) between a and b. Hence: 

fc-i 

dc,iP){a,b) =^\pjPj+i\. 
i=i 

Now, by proposition 14. 1[ for each edge {pj,pjj^i) there is a path in Si between pj and pj+i whose 
length is at most t ■ \pjPj+i\. Therefore: 

k-l k-1 

ds,{P)ia, b) <^t- \pjPj+i\ = t \pjPj+i\ = t ■ dc^^p){a, b) 

which means that in Si, there exists a path between a and b whose length is at most t ■ dci(p){a, b). 
□ 



Corollary 4.3 If S is a strong t-spanner of C (P), then S'nUDG(P) is a strong t-spanner of UT>G{P). 
Proof: Just notice that UDG = Ci for some i and the result follows from Proposition I4.2[ □ 

Thus, we have shown sufficient conditions for a graph to be a spanner of the unit disk graph. 
We now show that these conditions are also necessary. 

Proposition 4.4 If S is a subgraph of C{P) such that for all i = 1. . . (^), Si is a t-spanner of Ci{P), 
then S is a strong t-spanner of C{P). 

Proof: Let a, b be any pair of points chosen in P. We have to show that in S, there is a path 
between a and b such that 

1. its length is at most t ■ \ab\ and 
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2. every edge on the path has length at most \ab\. 

Let Ci = (a, b). We know that Si is a t-spanner of Ci(P). Since Ci{P) contains e^, dc^(p){a, b) = 
\ab\. Hence, there is a path in Si (and therefore in S) whose length is at most t ■ dc^[p){a, b) = t\ab\. 
Also, since it is in Si, all of its edges have length at most Li = \ab\. □ 




Figure 8: The ^-graph is not a strong t-spanner. 

The two last results, together, allow us to determine whether or not given families of geometric 
graphs are also spanners of the unit disk graph. First, since the graph is a strong t-spanner, we 
already know that it is also a spanner of the unit disk graph. Second, Bose et al. [ 2] showed that 
the Yao graph [7] and the Delaunay triangulation are strong t-spanners. Therefore, these graphs 
are also spanners of the unit disk graph. Third, the 0-graph[5] is not always a spanner of the unit 
disk graph. The reason for that is that in a cone, the edge you chose may not be the shortest edge. 
Hence, the path from a point p to a point q may contain edges whose length is greater than \pq\ 
(see Figure [H) . Using Proposition I4.4[ we thus know that the intersection of the ^-graph with the 
unit disk graph may not be a spanner of the unit disk graph. Indeed, the intersection of the 0-graph 
with the unit disk graph may not even be connected. 

5 Simulation Results 

In section EJ we provided worst-case analysis of the spanning ratio of the G\ graph. Using simu- 
lation, we now provide estimates of the average spanning ratio of the graph. Using a uniform 
distribution, we generated 200 sets of 200 points each and computed the spanning ratio for A rang- 
ing from 0.5 to 1 and 6 ranging from 5° to 90° (for 6 = 0°, the spanning ratio is exactly 1). For 
each graph, we then computed the spanning ratio and the local routing ratio. The spanning ratio is 
defined as the maximum, over all pair of points {p, q), of the length of the shortest from pto q path 
in the G^ graph divided by \pq\. The local routing ratio is defined as the maximum, over all pair of 
points {p, q), of the length of the path produced by using a local routing strategy in the G^ graph 
divided by \pq\. The local routing strategy we have used is the following: at each step, send the 
message to the neighbor which destroyed q. We also tried the strategy which consists in choosing 
the neighbor which is the nearest to q, and the results we obtained were the same. 

Figure [9] and Table [T] show the results we obtained for the spanning ratio. Figure [10] and 
Table [2] show the results we obtained for the local routing ratio. For the spanning ratio, the 95% 
confidence interval for these values is ±0.0319. For the local routing ratio, the 95% confidence 
interval is ±0.0735. 
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One interesting conclusion we can draw from these results is that for the spanning ratio, 9 has a 
more decisive influence than A. Figure [TTl shows the simulation results for the cases where A = 0.75. 
We see that even though both ratios generally increase when 9 increase, the spanning ratio varies 
between 1.07 and 2.21 (107% variation), while the local routing ratio only varies between 2.33 
and 2.77 (19% variation). For the local routing ratio, it is the other way around. It is A which has 
a more decisive influence. Figure [T2l shows the influence of A when 9 = 45°. In that case, the local 
routing ratio varies between 1.72 and 4.55 (165% variation), while the spanning ratio only varies 
between 1.52 and 1.81 (19% variation). 

6 Conclusion 

We conclude with the problem that initiated this research. Determine whether or not HSP is a 
strong t-spanner for some constant t. 
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Figure 9: Spanning Ratio for A = 0.5 to 1 and 6i = 5° to 90°. 
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1 


88 
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2.47 
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1.78 


1 
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2.10 


2.19 


2.27 


2.35 
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2.47 


2.51 
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1.78 


1 
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2.09 


2.21 


2.26 


2.33 


2.41 


2.46 


2.54 



Table 1: Spanning Ratio for A = 0.5 to 1 and = 5° to 90°. 
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Figure 10: Local Routing Ratio for A = 0.5 to 1 and 6* = 5° to 90°. 
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Table 2: Local Routing Ratio for A = 0.5 to 1 and = 


5° to 90°. 
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